Ground state of the parallel double quantum dot system 
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We resolve the controversy regarding the ground state ol the parallel double quantum dot system near hall 
filling. The numerical renormalization group (NRG) predicts an underscreened Kondo state with residual spin- 
1/2 magnetic moment, In 2 residual impurity entropy, and unitary conductance, while the Bethe Ansatz (BA) 
solution predicts a fully screened impurity, regular Fermi-liquid ground state, and zero conductance. We cal- 
culate the impurity entropy of the system as a function of the temperature using the hybridization-expansion 
continuous-time quantum Monte Carlo technique, which is a numerically exact stochastic method, and find ex- 
cellent agreement with the NRG results. We show that the origin of the unconventional behavior in this model 
is the odd-symmetry "dark state" on the dots. 

PACS numbers: 72.10.Fk, 72.15.Qm 



Quantum dot (QD) nanostructures serve as model systems 
for studying fundamental many-particle effects, such as the 
competition between the Kondo screening and the exchange 
interaction. These effects can give rise to ground states of 
the Fermi liquid or non-Fermi liquid nature. The fingerprints 
of those different states as well as the quantum-phase transi- 
tions between them have been predicted theoretically within 
the generalizations of the Anderson impurity model and ob- 
served in transport experiments H - 14 ] . 

To account for such a rich behavior powerful non- 
perturbative theoretical tools must be used. Among these 
techniques, the numerical renormalization group (NRG) 115 



Hit] is popular because of its wide applicability, reliability, 
and relatively low computational demands. The NRG results 
for the conductance of single QDs in the Kondo regime have 
played the key role in conclusively proving the occurrence of 
the Kondo effect in these systems yl [171 ll8[] . The NRG is an 
approximate method, expected to be asymptotically exact on 
the lowest energy scales. In constructing the effective chain 
Hamiltonian, the conductance band continuum is discretized, 
a decomposition into Fourier modes is performed in each in- 
terval, and a single representative state from each interval 
is retained. This approximation is, however, well controlled 
Il5tl . Some impurity models are integrable and can be solved 
exactly using the Bethe Ansatz (BA) technique; these ana- 
lytical solutions are very valuable as reference results which 
serve as benchmark for more generally applicable methods. 
The single-impurity Kondo and Anderson models are both in- 
tegrable | Ijl 2(1 and an excellent agreement has been found 
between the NRG and BA results 11911 . 

One of the simplest problems that exhibits non-Fermi liquid 
behavior is that of the two impurities coupled to conduction 
bands: the double quantum dot (DQD). Recently, a BA solu- 
tion has been proposed for a family of two-impurity models 
of DQDs HI- 24]. For two QDs coupled in parallel between 
the conductance leads with equal hybridization strengths, thus 
forming a symmetric ring (Fig. [T]), the problem maps onto a 



two-impurity model with a single effective conduction band. 
The BA solution predicts that near the particle-hole symmet- 
ric point the two electrons residing on the dots form a regular 
Fermi-liquid (FL) singlet state with the conduction band elec- 
trons, that the system has zero conductance at the particle-hole 
symmetric point, and that the standard Friedel sum rule is sat- 
isfied 12111 , All these predictions are, however, at odds with 
the NRG results JHQ. 



On general physical grounds one expects that at some high- 
energy scale the two impurity moments bind into a spin-triplet 
effective state due to the presence of the ferromagnetic RKKY 
interaction (the model corresponds to the r = limit of 
the standard two-impurity model J34-38]), then this spin is 
partially screened in the single-channel spin-1 Kondo effect 
yielding a singular Fermi liquid ground state, a residual spin- 
1/2 magnetic moment and residual In 2 entropy 1 3^- 41 ] . The 
NRG results are fully consistent with this scenario. This is 
also in line with the conventional wisdom that a single screen- 
ing channel can screen one half-unit of the impurity spin 
Furthermore, the system reaches unitary con- 



ductance at zero temperature. The Friedel sum rule in its stan- 
dard form 6 = (7r/2)ni mp is violated (there is an additional 
phase shift by 7r/2, see Ref.lill Sec. IV.B). The underscreened 
Kondo effect has already been experimentally observed in sys- 
tems described by two-orbital impurity models similar to the 
one discussed in this work 1 12lll3l 14211 . 



It is disconcerting that two purportedly highly reliable 
methods produce opposite results for an elementary impurity 
model. If the NRG were found to be flawed, this would de- 
mand a reinvestigation of the applicability of the method and 
put in question the reliability of a large number of published 
theoretical results. It has been suggested that the presumed 
deficiency of the NRG consists in disregarding the higher 
conduction-band modes in each discretization interval and a 
modified discretization has been proposed 12311 . Since the 
modified discretization scheme again maps onto a single ef- 
fective screening channel, the residual moment would still not 
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Figure 1 : Schematic representation of the model for a parallel double 
quantum dot coupled to two semi-infinite tight-binding chains with 
equal hopping constants. The Fermi level is fixed at zero energy. 



be screened, thus this does not solve the observed discrepancy. 

For this reason, in this work we resolve the contro- 
versy using an independent method: we perform extensive 
numerical simulations of the impurity Hamiltonian using 
the hybridisation-expansion continuous-time quantum Monte 
Carlo (CTQMC) algorithm JH-jH]. This method is numer- 
ically exact, its accuracy being limited solely by the calcula- 
tion time. No approximations nor simplifications of the model 
Hamiltonian need to be performed. The price for being exact 
are, however, heavy computational demands. We show that 
the simulated results are consistent with the NRG calculation, 
thus the proposed BA solution is not correct. 

The Hamiltonian is H = -ffband + ^dots + -Hhyb- Here 
-Hband = J2kaj £ kc\ aj c kaj is the conduction-band Hamilto- 
nian where k is momentum, a =t, 4- is spin, and j = 1,2 



indexes the leads. H, 



dots 



Un 



is the quan- 



tum dot Hamiltonian. The number operator n.j CT is defined as 
n ia = djadia and m — J2 a n i^> £ i s tne on-site energy, and 
U is the electron-electron repulsion. Finally, 
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is the coupling Hamiltonian, where L is a normalization con- 
stant. We model the conductance leads by semi-infinite tight- 
binding chains and the dots couple to the end of these chains 
by a hopping term, see Fig. Q] The hybridization function is 
defined as T(u) = 2np(uj)V 2 with p(ui) the density of states 
in the band and Vk the coupling coefficient at momentum k 
that corresponds to energy w, the additional factor 2 in the ex- 
pression takes into account that there are two leads. We have 
e = 2t cos k and Vk = f sin k, which gives 



r(w) = ry 1 - (w/2t) 2 , 



(2) 



with r = (t') 2 /t. In the following, we use the half-bandwidth 
D = 2t as the energy unit. The final expression for the impu- 
rity action is 



S= / dr 
Jo 



d ™ ( J~ - v + 6 ) di ° + u n ^ n ^ 

i=(l,2),cr 



- f dr f dr'dl(T)r(T-r')d e(T (r'). 
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The even/odd combination of operators is defined as d\, a = 

(d\ a ± d^J) /V%- Only the even orbital hybridizes with the 
conduction band. 

In this problem, there are three important energy scales 
1 2911 . On the scale of U the local moments are formed and 
the dots start to behave as two spin- 1/2 impurities. On the 
RKKY scale of Jrkky ~ U{pJ K ) 2 = {M/tt 2 )T 2 /U the 
spins bind into a S = 1 state. Here pJx = 8T/ttU quantifies 
the strength of the exchange coupling. Finally, on the scale 
of the Kondo temperature Tjc ~ Uy/pJx cxp (—1/pJx), the 
impurity moment is partially screened from spin-1 to spin-1/2 
in the single-channel spin-1 Kondo effect |39ll . No other low- 
energy scales are present in this problem [34]. The U /Y ratios 
in experiments on quantum dots range roughly from 1 to 20, 
thus the three energy scales introduced above are not neces- 
sarily well separated. 

Different methods for solving impurity models are best 
compared by calculating thermodynamic functions such as 
energy or entropy, since all other quantities of interest can 
be obtained by taking appropriate derivatives. In this work 
we calculate the impurity entropy Sj mp = S — , where 
S is the entropy of the full system, while S^°> is the entropy 
of the conduction band alone. Using the NRG one can com- 
pute Simp over a wide range of temperature scales with lit- 
tle numerical effort. To the contrary, the CTQMC becomes 
increasingly numerically demanding at low temperatures and 
for larger hybridization strength Y. The comparison of the re- 
sults can thus only be performed in a limited temperature win- 
dow which depends on Y. To circumvent this limitation, we 
perform calculations for a range of Y at fixed U; in this way 
we tune the characteristic low-energy scales of the problem 
(the Kondo temperature, the RKKY scale) from very small 
to very large values, making them pass through the available 
temperature window for different values of Y. 

In NRG, the impurity entropy is computed in the standard 
way iflrjll . To obtain good results even on the temperature 
scale of the bandwidth, the discretization is performed with 
a small value of the discretization parameter A = 1.8 and the 
twist-averaging with N z = 16 is used. Nevertheless, some 
small quantitative systematic errors due to the discretization of 
the band are expected for temperatures approaching the band- 
width. The truncation is performed with a sufficiently high 
energy cutoff E cuto g = 12lon that the results can be consid- 
ered as fully converged. 

In CTQMC, we calculate the impurity energy for a range of 
temperatures and obtain the entropy as 
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The integration constant is fixed by the high-temperature 
asymptotic limit of 2 In 4. The impurity energy is defined as 
E c with 
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The notation ()o indicates that the expectation value is com- 
puted for the system without the impurity (i.e., in the r — > 
limit). E c is obtained from the impurity spectral function in 
the Matsubara space Q(iui n ): 



E c = — d(iuj n )g(iuj n ), 



where 



d(z) 



-isgn[Im(z)] 



(6) 



(7) 



The impurity spectral function is defined as Q = Qn + Qyi + 
Qi\ + Q22 with Qij(z) = {{di;dj)) z . To avoid the minus 
sign problem, the simulation is performed in the even/odd ba- 
sis. Since the asymptotic behavior of the Green's function is 
Q(iu> n ) ~ J— , we subtract and add l/(iuj n ) and calculate the 
problematic part exactly. 
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Figure 2: (Color online) Comparison of the impurity entropy curves 
at the particle-hole symmetric point. Lines: NRG, circles: QMC 
(with error bars). Energy unit is the half-bandwidth D = 1. 
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Figure 3: (Color online) Comparison of the impurity entropy curves 
away from the particle-hole symmetric point. Lines: NRG, circles: 
QMC (with error bars). 

The results for the temperature dependence of the impu- 
rity entropy are shown in Fig. [2] at the p-h symmetric point 



and in Fig.[3]away from it. The agreement between the NRG 
and CTQMC results is excellent in all cases considered. A 
small systematic deviation of the NRG from the QMC results 
at very high temperatures (T ~ D) is anticipated. At interme- 
diate temperatures, the agreement improves until at some low 
T-dependent temperature the QMC simulation can no longer 
be performed in reasonable time due to slow thermalization. 
Nevertheless, NRG and QMC are found to agree below all 
the relevant energy scales in the problem. The DQD system 
near the p-h symmetric point thus behaves as a singular FL, as 
predicted by the NRG. Furthermore, in Fig.[3]we show numer- 
ical evidence of a quantum phase transition (QPT) where, as a 
function of e, the system goes from a singular FL to a regular 
FL ground state (the BA solution predicts no such transition). 




Figure 4: (Color online) Phase diagram in the (e, U) plane obtained 
using the NRG. The entropy curve at point A is shown in Fig. [2] 
while those at points B and C are compared in Fig. [3] 



The DQD problem being integrable, this raises the ques- 
tion why the BA solution differs from the NRG and QMC 
results. One possibility is that the BA wave-function cor- 
responds to some excited state instead of the actual ground 
state of the system. We remark that the NRG and QMC are 
both grand-canonical-ensemble calculations in the thermody- 
namic limit, thus the occupancy of the dots is automatically 
correctly determined, while in BA the thermodynamic limit 
is taken at the end of the calculation, thus one needs to take 
care to choose the wavefunction from the correct charge and 
spin sector. In particular, the occupancy of the odd state d\ a 
is important in this problem. In the non-interacting limit, one 
has [if, n odd ] = with n odd = J^a dt. a d . a , thus the oc- 
cupancy of the odd state is a conserved quantity (the system 
is non-ergodic). This state is completely decoupled from the 
continuum and is sometimes referred to as the bound state in 
the continuum, dark state, or ghost Fano resonance |3(| 46Tl . 
It lies exactly at the Fermi level when the system is tuned to 
the particle-hole symmetric point. Formally, at this point the 
system has 2 In 2 residual entropy since for each spin the level 
may be either occupied or unoccupied at no energy cost. As 
the interaction is turned on, n dd is no longer a constant of 
motion, yet the odd state still plays a non-trivial role. We find 
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that the "dark state" gives rise, at finite U, to a finite range of 
the gate voltages e around the particle-hole symmetric point 
where the occupancy of the odd level is pinned exactly to one, 
i.e., the odd level accommodates an electron of either spin, 
which is asymptotically decoupled from the rest of the sys- 
tem on low energy scales. This gives rise to a singular FL 
ground state with In 2 residual entropy. The breakdown of 
the regular FL ground state in this parameter range has been 
noted previously using the Gunnarson-Schonhammer varia- 
tional approach based on a regular FL trial function, which 
failed to converge in the parameter range that is now known 
to correspond to the singular FL phase 047I1 . The correspond- 
ing interval of the gate voltages e is delimited by QPTs across 
which the impurity charge changes discontinuously (in both 
even and odd levels) |29[ 3(1. These QPTs emerge directly 
from the non-interacting dark state as the interaction is turned 
on, see Fig. [4] We also note that the system is fully conduct- 
ing at the particle-hole symmetric point for any value of U, 
including U = 0. It remains an open technical problem how 
to include these effects in the BA calculation in order to obtain 
the correct doublet ground state around half filling. 

In conclusion, we solved the problem of the parallel DQD 
using hybridization expansion CTQMC and demonstrated that 
the results agree with the results of the NRG. We pointed out 
the non-trivial role played by the "dark state" odd orbital. We 
remark that the studied model is an excellent benchmark test 
for various many-body techniques, since only few appear to 
produce the correct results, and that even supposedly exact an- 
alytical calculations need to be a-posteriori validated against 
other reference results. 
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